The Monte Carlo Event Generator DPMJET-III* 



Stefan Roesler 1 , Ralph Engel 2 , and Johannes Ranft 3 

1 SLAC, P.O. Box 4349, Stanford CA 94309, USA 

2 University of Delaware, Bartol Res. Inst., Newark DE 1 

3 University of Siegen, D-57068 Siegen, Germany 



Abstract. A new version of the Monte Carlo event generator Dpmjet is presented. 
It is a code system based on the Dual Parton Model and unifies all features of the 
Dtunuc-2, Dpmjet-II and Phojet1.12 event generators. Dpmjet-III allows the sim- 
ulation of hadron-hadron, hadron-nucleus, nucleus-nucleus, photon-hadron, photon- 
photon and photon-nucleus interactions from a few GeV up to the highest cosmic ray 
energies. 



1 Introduction 

Hadronic collisions at high energies involve the production of particles with low 
transverse momenta, the so-called soft multiparticle production. The theoretical 
tools available at present are not sufficient to understand this feature from QCD 
and phenomenological models are typically applied instead. The Dual Parton 
Model (DPM) jl| is such a model and its fundamental ideas are presently the 
basis of many of the Monte Carlo (MC) implementations of soft interactions in 
codes used for Radiation Physics simulations. 

Many of these implementations are however limited in their application by, 
for example, the collision energy range which they are able to describe or by 
the collision partners (hadrons, nuclei, photons) which the model can be used 
for. With respect to modern multi-purpose codes for particle interaction and 
transport these limitations at high energy are clearly often a disadvantage. 

In this paper we present the Dpmjet-III code system, a MC event generator 
based on the DPM which is unique in its wide range of application. Dpmjet-III is 
capable of simulating hadron-hadron, hadron-nucleus, nucleus-nucleus, photon- 
hadron, photon-photon and photon-nucleus interactions from a few GeV up to 
the highest cosmic ray energies. 

In the present paper we give an overview over the different components and 
models of Dpmjet-III and present a few examples for comparisons of model 
results with experimental data. 



2 The Concept of the Program 

Dpmjet-III is the result of merging all features of the event generators Dpmjet- 
II §U and Dtunuc-2 |,§ into one single code system. The latter two codes 
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are similar in their underlying concepts, however they differ in the Monte Carlo 
realization of these concepts, in particular, of the DPM. 

Whereas individual nuclcon-nucleon collisions in Dpmjet-II are simulated 
based on the Dtujet model B, Dtunuc-2 is using Phojet1.12 Since 
Phojet describes not only hadron-hadron interactions but also hadronic inter- 
actions involving photons, Dtunuc-2 allows also the simulation of photopro- 
duction off nuclei. Therefore, the strength of Dtunuc-2 is in the description of 
photoproduction and nuclear collisions up to TeV-energies. On the other hand, 
Dpmjet-II is widely used to simulate cosmic-ray interactions up to the highest 
observed energies ||. 

However, many program modules in Dpmjet-II and Dtunuc-2 are also iden- 
tical. Examples are the Glauber-Gribov formalism for the calculation of nuclear 
cross sections |0L the formation-zone intranuclear cascade , the treatment of 
excited nuclei and the HADRiN-model for the description of interactions 

below 5 GeV |f] . 

The core of Dpmjet-III consists of Dtunuc-2 and Phojet1.12. In addition 
all those features of Dpmjet-II were added which were not part of Dtunuc-2 
so far. This includes, for example, quasi-elastic neutrino interactions |l4| and 
certain baryon-stopping diagrams |l5t . 

3 Models Implemented in DPMJET-III 

3.1 The Realization of the Dual Parton Model 

The DPM combines predictions of the large N c , Nf expansion of QCD |l6) and 
assumptions of duality with Gribov's reggeon field theory Phojet, 
being used for the simulation of elementary hadron-hadron, photon-hadron and 
photon-photon interactions with energies greater than 5 GeV, implements the 
DPM as a two-component model using Reggeon theory for soft and perturba- 
tive QCD for hard interactions. In addition to the model features as described in 
detail in the version 1.12 incorporates a model for high- mass diffraction dis- 
sociation including multiple jet production and recursive insertions of enhanced 
pomeron graphs (triple-, loop- and double-pomeron graphs). In the following 
only the new features are briefly discussed. 

High-mass diffraction dissociation is simulated as pomeron-hadron or pomeron- 
pomeron scattering, including multiple soft and hard interactions [p0[ . To ac- 
count for the nature of the pomeron being a quasi-particle, the CKMT pomeron 
structure function |2l| with a hard gluonic component is used. These consid- 
erations refer to pomeron exchange reactions with small pomeron-momentum 
transfer, \t 2 \. For large \t 2 \ the rapidity gap production (e.g. jet-gap-jet events) 
is implemented on the basis the color evaporation model |2^] . 

Extrapolating the two-channel eikonal-unitarization of a hadron-hadron am- 
plitude as used in Phojet to very high energies raises the question of the treat- 
ment of enhanced graphs which become more and more important at high energy 
and lead to large multiplicity fluctuations. A full amplitude calculation including 
enhanced graphs is very involved and not suited for a Monte Carlo implemen- 
tation. Therefore, based on the results of |23], we use the simpler approach of 
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interpreting each soft pomeron as the sum of a series of a bare soft pomeron 
and enhanced graphs (Froissaron). In practice, this results in the simulation 
of possibly recursive subdivisions of a single Froissaron cut into various other 
configurations such as, for example, two cut pomerons or a single cut pomeron 
and a diffractive scattering. However, the current implementation should only 
be considered as a first step toward a consistent treatment of enhanced graphs 
at very high energy because of its limitation to soft interactions. 

3.2 Hadronic Interactions Involving Photons 

The photon is assumed to be a superposition of a bare photon interacting in 
direct processes and a hadronic photon interacting in resolved processes. 

The description of interactions of the hadronic photon with nuclei is based 
on the Generalized Vector Dominance Model (GVDM) p4| . Photons are as- 
sumed to fluctuate into quark-antiquark states V of a certain mass M and the 
interaction is described as scattering of the hadronic fluctuation on the nucleus. 
Correspondingly, the scattering amplitude ay a reads Q 



a VA (s,Q 2 ,M 2 ,B) = / [[al 3 ^ i>* A a VA (s,Q 2 ,M 2 ,B 1 ,...,B A ) ^ A (1) 



where ay a is expressed in terms of interactions on individual nucleons N ac- 
cording to the Gribov-Glauber picture (see below). The model is limited to low 
photon-virtualities Q 2 satisfying the relation Q 2 <C 2rriNV (v and mjv being the 
photon energy and nucleon mass). For individual gg- nucleon interactions it is 
sufficient to consider only two generic qq-st&tes, the first one grouping p°, to and 
4> and 7r + 7r _ -states up to the 0-mass together and the second one including all 
gg-states with higher masses 0. 

Direct photon interactions are treated as either gluon-Compton scattering or 
photon-gluon fusion processes on a single nucleon. The consideration of so-called 
anomalous interactions allows a steady transition between direct and resolved 
interactions [Q. 

Finally, an interface to Lepto6.5 |^5) allows to simulate deep-inelastic scat- 
tering off nuclei. 

3.3 The Gribov-Glauber Multiple Scattering Formalism 

The Monte Carlo realization of the Gribov-Glauber multiple scattering formal- 
ism follows the algorithms of || and allows the calculation of total, elastic, 
quasi-elastic and production cross sections for any high-energy nuclear collision. 
Parameters entering the hadron-nuclcon scattering amplitude (total cross section 
and slope) are calculated within Phojet. 

For photon-projectiles ideas of the GVDM have been incorporated in order 
to correctly treat the mass of the hadronic fluctuation and its coherence length 




avA{s 1 Q 2 ,M 2 ,B 1 



>Ba) = \ ( 1 - A [1 + 2ia VN (s, Q 2 ,M 2 , B v j\ I (2) 
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as well as pointlike photon interactions Q. Realistic nuclear densities and radii 
are used for light nuclei and Woods-Saxon densities otherwise. 

During the simulation of an inelastic collision the above formalism sam- 
ples the number of "wounded" nucleons, the impact parameter of the colli- 
sion and the interaction configurations of the wounded nucleons. Individual 
hadron(photon,nucleon)-nucleon interactions are then described by Phojet in- 
cluding multiple hard and soft pomeron exchanges, initial and final state radia- 
tion as well as diffraction. 

As a new feature, Dpmjet-III allows the simulation of enhanced graph cuts 
in non-diffractive inelastic hadron-nucleus and nucleus-nucleus interactions. For 
example, in an event with two wounded nucleons, the first nucleon might take 
part in a non-diffractive interaction whereas the second one scatters diffracively 
producing only very few secondaries. Such graphs are prediced by the Gribov- 
Glauber theory of nuclear scattering but are usually neglected. 

Finally, all color neutral strings are hadronizcd according to the Lund model 
as implemented in Pythia f|(||27j . 



3.4 The Intranuclear Cascade and Break-up of Excited Nuclei 

The treatment of intranuclear cascades in spectator prefragments and their sub- 



sequent fragmentation is largely identical to the one described in Refs. 1 11 , 12j] . 

Particles created in string fragmentation processes are followed on straight 
trajectories in space and time. A certain formation time is required before newly 
created particles can re-interact in the spectator nuclei. These re-interactions are 
of low energy and are described by Hadrin[|i3| based on parameterized exclusive 
interaction channels. In nucleus- nucleus collisions the intranuclear cascade is 
calculated in both the projectile and target spectators. 

Excitation energies of prefragments are calculated by summing up the recoil 
momenta transfered to the respective prefragment by the hadrons leaving the 
nuclear potential (a constant average potential is assumed). The prefragments 
are assumed to be in an equilibrium state and excitation energy is dissipated by 
the evaporation of nucleons and light nuclei and by the emission of photons. 



4 Comparison to Experimental Data 

Since Dpmjet-III is the result of merging Dpmjet-II and Dtunuc-2 its predic- 
tions have to be in agreement to experimental data where there was agreement 
for the two latter codes before. However, this has to be proven again. Here, only a 
few examples are given which should represent the large amount of comparisons 
of Dpmjet-III results with experimental data which exist. 

Fig. |l|a shows the transverse momentum distribution of negative hadrons from 
p-W collisions together with data (28). The rapidity distributions of negative 
hadrons in central S-S and S-Ag collisions are compared to data | f29[ in Fig. |l|b. 
Two examples for interactions involving photons are given in Fig. |2j. Hadronic 
interactions of muons are described by the radiation off the muon of a quasi- 
real photon and the subsequent interaction of the photon. Fig. 0a shows average 
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Fig. 1. Negatively charged hadron production in nuclear collisions at 200 GeV/nucleon 

multiplicities of charged hadrons from /x-Xe interactions at 490 GeV compared to 
data JhJ. In Fig. [2}d the calculated inclusive transverse momentum cross section 
of charged particles produced in two-photon collisions at LEP is compared to 
the combined data set of the ALEPH, L3, and OPAL Collaborations for low-Q 2 
deep inelastic scattering |pl[. 
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Fig. 2. Comparison of Dpmjet-III results to data on interactions involving photons. 

5 Conclusions 

A new version of the Dpmjet event generator is presented. Dpmjet-III is based 
on Dpmjet-II, Dtunuc-2 and Phojet1.12 and unifies all features of these three 
event generators in one single code system. It has been thoroughly tested and, 
in due time, will largely superseed the older Dpmjet and Dtunuc versions. 

It is presently not advisable to use the code for very low-energy nucleus- 
nucleus collisions (below w 10—20 GeV). This requires further testing and tuning 
of parameters. Furthermore deficiencies exist in the description of some effects 
observed in heavy ion collisions at AGS- and SPS-energies (e.g. strangeness 
enhancement, transverse energy flow). 

The code is available on request from the authors (Stefan. Roesler@cern.ch, 
Johannes. Ranft@cern.ch) and further information can be found on the World 
Wide Web (http://home.cern.ch/sroesler/dpmjet3.html). 



6 Stefan Roesler et al. 



6 Acknowledgements 

The work of S.R. and R.E. is supported by the Department of Energy under 
contracts DE-AC03-76SF00515 and DE-FG02-91ER40626, respectively. 

References 

1. A. Capella et al.: Phys. Rept. 236, 225 (1994) 

2. J. Ranft: Phys. Rev. D 51, 64 (1995) 



3. J. Ranft: Dpmjet version II. 5: Code manual, hep-ph/9911232 , Preprint University 
of Siegen SI-99-6 (1999) 

4. R. Engel, J. Ranft and S. Roesler: Phys. Rev. D 55, 6957 (1997) 

5. S. Roesler, R. Engel and J. Ranft: Phys. Rev. D 57, 2889 (1998) 

6. P. Aurenche et al.: Phys. Rev. D 45, 92 (1992) 

7. R. Engel: Z. Phys. C 66, 203 (1995) 

8. R. Engel and J. Ranft: Phys. Rev. D 54, 4244 (1996) 

9. S.Y. Shmakov, V.V. Uzhinskii and A.M. Zadorozhny: Comput. Phys. Commun. 54, 
125 (1989) 

10. J. Ranft: Phys. Rev. D 37, 1842 (1988) 

11. A. Ferrari, J. Ranft, S. Roesler and P.R. Sala: Z. Phys. C 70, 413 (1996) 

12. A. Ferrari, J. Ranft, S. Roesler and P.R. Sala: Z. Phys. C 71, 75 (1996) 

13. K. Hanfigen and J. Ranft: Comput. Phys. Commun. 39, 37 (1986) 

14. G. Battistoni, P. Lipari, J. Ranft and E. Scapparone: Simulation of nuclear effects 



in quasielastic and resonant neutrino interactions, hep-ph/9801426, Preprint INFN- 
AE-03-98 (1998) 

15. J. Ranft, R. Engel and S. Roesler: Baryon stopping in high energy collisions in the 
Dpmjet-III model, Presented at the same meeting by J. Ranft 

16. G. Veneziano: Nucl. Phys. B 74, 365 (1974) 

17. G. F. Chew and C. Rosenzweig: Phys. Rept. 41, 263 (1978) 

18. V. N. Gribov: Sov. Phys. JETP 26, 414 (1968) 

19. R. Engel: H adronic interactions of photons at h igh energies, Ph.D. thesis, Univer- 



sitat Siegen, ittp : //lepton.bartol .udel . edu/~eng/phojet .html (1997) 
20. F. W. Bopp, R. Engel and J. Ranft: Rapidity gaps and the PHOJET Monte Carlo, 



^iep-ph/9803437 , in Proceedings of LAFEX Int. School on High-Energy Physics 
(LISHEP98), Session C: Workshop on Diffractive Physics, Rio de Janeiro, Brazil, 
February 16-20, 1998 

21. A. Capella et al.: Phys. Rev. D 53, 2309 (1996) 

22. O. J. P. Eboli, E. M. Gregores and F. Halzen: Phys. Rev. D 58, 114005 (1998) 

23. A. B. Kaidalov, L. A. Ponomarev and K. A. Ter-Martirosyan: Sov. J. Nucl. Phys. 
44, 468 (1986) 

24. A. Donnachie and G. Shaw: Generalized Vector Dominance, in Electromagnetic 
Interactions of Hadrons, Volume 2, ed. by A. Donnachie and G. Shaw (Plenum 
Press, New York 1978) 

25. G. Ingelman, A. Edin and J. Rathsman: Comput. Phys. Commun. 101, 108 (1997) 

26. T. Sjostrand: Comput. Phys. Commun. 82, 74 (1994) 



27. T. Sjostrand: Recent progress in PYTHIA, |hep-ph/0001032j , Preprint Lund Uni- 
versity LU TP 99-42 (1999) 

28. The HELIOS Collaboration: T. Akesson et al.: Z. Phys. C 46, 361 (1990) 

29. The NA35 Collaboration: T. Alber et al.: Eur. Phys. J. 2, 643 (1998) 

30. The E665 Collaboration: M. R. Adams et al.: Z. Phys. C 61, 179 (1994) 

31. A. J. Finch: Nucl. Phys. B (Proc. Suppl.) 82, 156 (2000) 



